a^ 



X) 



0^ 



1. INTRODUCTION 



Cosmic antiprotons as a probe 
for supersymmetric dark matter? 

Lars Bergstrom*, Joakim Edsjo^, Piero Ullio"'' 

Department of Physics, Stockholm University, Box 6730, SE-113 85 Stockholm, Sweden 

February 1, 1999 

ABSTRACT 

The flux of cosmic ray antiprotons from neutralino annihilations in the galactic halo is computed for 
a large sample of models in the MSSM (the Minimal Supersymmetric extension of the Standard Model). 
We also revisit the problem of estimating the background of low-energy cosmic ray induced secondary 
antiprotons, taking into account their subsequent interactions (and energy loss) and the presence of 
nuclei in the interstellar matter. 

We consider a two-zone diffusion model, with and without a galactic wind. We find that, given 
Q«J ' the uncertainties in the background predictions, there is no need for a primary (exotic) component to 

O^ . explain present data. However, allowing for a signal by playing with the uncertainties in the background 

estimate, we discuss the characteristic features of the supersymmetric models which give a satisfactory 
description of the data. We point out that in some cases the optimal kinetic energy to search for a signal 
^ . from supersymmetric dark matter is above several GeV, rather than the traditional sub-GeV region. 

Ft , ' The large astrophysical uncertainties involved do not, one the other hand, allow the exclusion of any 

of the MSSM models we consider, on the basis of data. 

We present besides numerical results also convenient parameterizations of the antiproton yields of all 
'basic' two-body final states. We also give examples of the yield and differential energy spectrum for a 
►^ , set of supersymmetric models with high rates. 

fvj ' We also remark that it is difficult to put a limit on the antiproton lifetime from present measurements, 

since the injection of antiprotons from neutralino annihilation can compensate the loss from decay. 

o 

^P , The mystery of the dark matter in the Universe remains unsolved. Among the most plausible candidates are Weakly 
Interacting Massive Particles (WIMPs), of which the supersymmetric neutralino is a favourite candidate from the point of 
view of particle physics. The neutralino arises naturally in supersymmetric extensions of the standard model, and has the 
attractive feature of giving a relic density which in large regions of parameter space is adequate to explain cosmological 
Q^' dark matter. We will in this paper consider the neutralino as a dark matter candidate within the Minimal Supersymmetric 
extension of the Standard Model (MSSM). For a thorough review of neutralino dark matter, see Jungman, Kamionkowski 
y^ • & Griest (1996). 

"^ ' Neutralino dark matter can be and is searched for in several ways: directly through detection of nuclear recoils and/or 
^ ionization in direct detection experiments, and indirectly through searches for their annihilation products from annihilation 
ILJ in the Earth or Sun (for neutrinos) and the galactic halo. In this paper, we discuss the detection prospects of antiprotons 
. j—{ from neutralino annihilation in the galactic halo. 

^^ As antimatter seems not to exist in large quantities in the observable Universe, including our own Galaxy, any contribu- 

$H ' tion to the cosmic ray generated antimatter flux (besides antiprotons also positrons) from exotic sources may in principle 
5t 1 be a good signature for such sources. Since neutralinos are constrained by supersymmetry to be Majorana fermions 
they are their own antiparticles and therefore the final state in their annihilations in the halo will contain equal amounts 
of matter and antimatter (given the particle physics constraints on CP violating couplings and in particular on baryon 
number violation). The excess of particles would drown in the background of particles from astrophysical sources, but 
there is a chance that antiparticles from this new primary source could be detectable. This issue has recently come into 
new focus thanks to upcoming space experiments like Pamela (Adriani & al. 1995) and Ams (Ahlen & al. 1994) with 
increased sensitivity to the cosmic antimatter flux. 

Cosmic ray induced secondary antiprotons are generated mainly through pp — > p + A" collisions of cosmic ray protons 
with interstellar matter. For kinematical reasons they are born with a non-zero momentum. The strategy to search for 
exotic signals has thus been to investigate the low-energy region since, e.g., a neutralino- induced component does not drop 
as fast at low energies. However, as we will see, this ideal picture is blurred to a large extent by a 'tertiary' component 
caused by scattering with energy loss of the secondary antiprotons. Also, heavier nuclei in the interstellar medium target 
(primarily helium) cause a significant antiproton flux at low energy. In particular, it is known that in proton-nucleus 
collisions antiprotons may be produced well below the nominal pp energy threshold. In addition, low energy particles have 
difficulties entering the heliosphere which makes the connection of the measured fluxes to the insterstellar ones dependent 
on a not completely known correction due to this solar modulation (which follows the 11-year solar cycle). 
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Using reasonable parameters for the computation of all these effects, we will show that we are able to explain in a 
satisfactory way the present experimental data on cosmic ray antiprotons without the need of a primary component. This 
is seen in Fig. 1, where our computed reference distribution for the background is compared to the recent data from the 
Bess experiment (Bess 95: Matsunaga & al. 1998; Bess 97: Orito 1998). The satisfactory agreement can be compared 
to the conclusions of a recent analysis (Bottino & al. 1998), where the need for an exotic component was more apparent. 
The main cause of this difference lies in our improved treatment of energy loss during propagation and our inclusion of 
the helium component of the interstellar medium. From Fig. 1 one can also see that the more recent Bess data (Orito 
1998) indicate a lower flux at low energy than earlier data. If reacceleration effects were included, the need for an exotic 
component would be even less. 

It is, however, also evident from Fig. 1 that the statistical sample of antiprotons presently available is very limited, 
so that a new primary component can not yet be ruled out with high significance even if the propagation parameters 
were known. By varying all parameters in this complex astrophysical problem, the room for an exotic contribution can, 
as we will see, be made quite large, so it is certainly worthwhile to investigate what a favoured dark matter candidate 
like the MSSM neutralino may yield in terms of a signal. It is apparent from Fig. 1 that the low-energy tail of cosmic 
ray induced antiprotons does not fall very rapidly with decreasing kinetic energy, so this particular region of phase-space 
may not yield as a nice signature for a dark matter neutralino as previously thought. Therefore, we will also discuss 
the cases (mainly for high-mass neutralinos) when the optimal kinetic energy for finding an exotic contribution is above 
several GeV. Unfortunately, it will turn out that the rates in that energy region are not large enough to cause a spectral 
distortion, unless the background of secondary antiprotons is considerably smaller and/or the signal is considerably larger 
than our canonical results show. This lack of spectral features in the neutralino-induced antiproton flux causes severe 
fundamental limitations for this indirect method of detecting supersymmetric dark matter. 

The idea of a dark matter-induced component in the cosmic antiproton spectrum has a long history. An early report 
(Buffington & al. 1981) of anomalous excess of cosmic ray antiprotons at low energies led to the suggestion that annihilation 
of relic neutralinos could be the source. Calculations of fluxes have since then been performed with different degrees of 
sophistication, ranging from order of magnitude estimates (Silk & Srednicki 1984), analytical expressions (Stecker, Rudaz 
& Walsh 1985) to results from Monte-Carlo simulations (Ellis & al. 1988; Stecker & Tylka 1989; Jungman & Kamionkowski 
1994; Bottino & al. 1995; Bottino & al. 1997). The latter method has also been used together with p production from 
neutralino annihilation in a minimal supergravity scheme (Diehl & al. 1995). 

In this work, we use the Lund Monte Carlo Pythia 6.115 (Sjostrand 1994) to simulate the energy spectrum of an- 
tiprotons from neutralino annihilation. We have used large numerical tables for our computations, but for convenience 
we also present useful parametrisations of the p fluxes from various annihilation channels. In the following Sections, we 
will describe the MSSM model we use, describe the Monte Carlo simulations, discuss the antiproton propagation model, 
discuss the background fluxes and the uncertainties in both the background and the signal. Finally, we will show and 
discuss our results. 
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Fig. 1. — The background of antiprotons solar modulated with (ftp ~ 500 MV. The Bess 95 and 97 data are also shown 
(Matsunaga & al. 1998; Orito 1998). 



Parameter 
Unit 


GeV 


M2 
GeV 


tan/3 
1 


niA 
GeV 


mo 
GeV 


Ab/ruQ 
1 


At/mo 
1 


Min 

Max 


-50000 
50000 


-50000 
50000 


1.0 
60.0 



10000 


100 
30000 


-3 
3 


-3 
3 



Table 1 

The ranges of parameter values used in our scans of the MSSM parameter space. Note that several special scans aimed at 
interesting regions of the parameter space has been performed. In total we have generated approximately 116 000 models 
which obey all accelerator constraints. Of these, about 41000 have a relic density in the range 0.025 < fty,h^ < 1. 

2. DEFINITION OF THE SUPERSYMMETRIC MODEL 

We work in the Minimal Supersymmetric Standard Model (MSSM). In general, the MSSM has many free parameters, 
but with some reasonable assumptions we can reduce the number of parameters to the Higgsino mass parameter /i, the 
gaugino mass parameter M2 , the ratio of the Higgs vacuum expectation values tan /3, the mass of the CP-odd Higgs boson 
mA (or mjjo), the scalar mass parameter ttiq and the trilinear soft SUSY-breaking parameters Ab and At for the third 
generation. In particular, we do not impose any restrictions from supergravity other than gaugino mass unification, which 
relates the other gaugino mass parameter Mi to M2. (We remind that one of the most attractive features of the MSSM is 
that, unlike the non-supersymmetric Standard Model, it is compatible with gauge coupling unification given the current 
data on the running of low-energy gauge couplings.) For a more detailed definition of the parameters and a full set of 
Feynman rules we refer to (Edsjo & Gondolo 1997; Edsjo 1997). 

The lightest stable supersymmetric particle is in most models the lightest neutralino (which we will henceforth just call 
'the neutrahno', x), which is a superposition of the superpartners of the gauge and Higgs fields, 

X = X? = NiiB + NuW"" + NuH^ + N^HI (1) 

It is convenient to define the gaugino fraction of the neutralino, 

Zg = \Nii\'' + \Ni2?- (2) 

For the masses of the neutralinos and charginos we use the one-loop corrections from the literature (Drees & al. 1997; 
Pierce & Papadopoulos 1994; Pierce & Papadopoulos 1994a; Lahanas, Tamvakis & Tracas 1994) and for the Higgs boson 
masses we use the leading logarithmic two-loop radiative corrections, calculated within the effective potential approach 
given by Carena & al. (1995). 

We make extensive scans of the model parameter space, some general and some specialized to interesting regions. 
In total we make 22 different scans of the parameter space. The scans are done randomly and are mostly distributed 
logarithmically in the mass parameters and in tan/3. For some scans the logarithmic scan in /i and M2 has been replaced 
by a logarithmic scan in the more physical parameters m~^ and Zg/ [1 — Zg) where m^ is the neutralino mass. Combining 
all the scans, the overall limits of the seven MSSM parameters we use are given in Table 1. 

We check each model to see if it is excluded by the most recent accelerator constraints, of which the most important 
ones are the LEP bounds (Carr 1998) on the lightest chargino mass, 

r 91 GeV , |m^+ - m^o| > 4 GeV 
^^t ^ \ 85 GeV , otherwise ' ^' 

and on the lightest Higgs boson mass tti^o (which range from 72.2-88.0 GeV depending on sin(/3 — a) with a being a 
mixing angle in the Higgs sector) and the constraints from b ^> sj (Ammar et al. 1993; Alam et al. 1995). The new 
higher-precision measurement from CLEG (Glenn & al. 1998) gives a slightly smaller range for that process than the one 
we have allowed; we have checked, however, that this causes no major changes in the properties related to p yield for the 
allowed models. 

For each allowed model we compute the relic density of neutralinos Vl~^h^ , where Q,^ is the density in units of the critical 
density and h is the present Hubble constant in units of 100 km s~^ Mpc^^. We use the formalism of Gondolo & Gelmini 
(1991) for resonant annihilations, threshold effects, and finite widths of unstable particles and we include all two-body 
tree- level annihilation channels of neutralinos. We also include the so-called coannihilation processes according to the 
results of Edsjo & Gondolo (1997) in the relic density calculation. 

Present observations favor h — 0.6 ± 0.1, and a total matter density VLm — 0.3 ± 0.1, of which baryons may contribute 
0.02 to 0.08 (see, e.g., Schramm & Turner 1998). Not to be overly restrictive, we accept Q,^h? in the range from 0.025 to 
1 as cosmologically interesting. The lower bound is somewhat arbitrary as there may be several different components of 
non-baryonic dark matter, but we demand that neutralinos are at least as abundant as required to make up the dark halos 
of galaxies. In principle, neutralinos with fl^h^ < 0.025 would still be relic particles, but only making up a small fraction 
of the dark matter of the Universe. We will consider models with Jl^ft,^ < 0.025 only when discussing the dependence of 
the signal on Vt^K^ . 
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021 
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032 


041 
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cc 


1.70 


1.40 










3.12 


0.04 


-2.22 





-0.39 


-0.076 








bb 


1.75 


1.40 










1.54 


0.11 


-2.22 





-0.31 


-0.052 








tt 


1.35 


1.45 










1.18 


0.15 


-2.22 





-0.21 











W+W~ 


306.0 


0.28 


7.2 


X 10-* 


2.25 


2.32 


0.05 


-8.5 


-0.31 


-0.39 


-0.17 


-2.0 X 10^2 


0.23 


Z''Z° 


480.0 


0.26 


9.6 


X 10"* 


2.27 


2.17 


0.05 


-8.5 


-0.31 


-0.33 


-0.075 


-1.5 X 10"" 


0.71 


99 


2.33 


1.49 










3.85 


0.06 


-2.17 





-0.312 


-0.053 









Table 2 

The parameters fitted to the antiproton distributions for neutrahno masses 50-5000 GeV. The parameters a23, a24, 033 
and 034 are always zero and not given in the table. The parameterized distributions are given by Eqs. (4)-(5). 

It may also be of interest to consider specifically models which naturally give a value of fl^h'^ which is close to the 
present 'best fit' value. We will therefore present some Figures where the range between 0.1 and 0.2 for fl^h^ is shown 
by special symbols. Since in general a small relic density implies a large annihilation cross section and vice versa, this 
tends to cut out a large fraction of the models with otherwise observable rates - a fact not always highlighted in previous 
analyses. 

3. ANTIPROTON PRODUCTION BY NEUTRALINO ANNIHILATION 

3.1. Introduction 

Neutralinos are Majorana fermions and will annihilate with each other in the halo producing Icptons, quarks, gluons, 
gauge bosons and Higgs bosons. The quarks, gauge bosons and Higgs bosons will decay and/or form jets that will give 
rise to antiprotons (and antineutrons which decay shortly to antiprotons). 

At tree level the relevant final states for p production will be qq, U, W+W^ , Z°Z°, W+H", ZH^, ZH^, H^H^ and 
H^H^- We will include all the heavier quarks (c, h and t), gauge bosons and Higgs boson final states in our analysis. In 
addition, we will include the Z7 (UUio & Bergstrom 1998) and the 2 gluon (Drees & al. 1994; Bergstrom & Ullio 1997) 
final states which occur at one loop-level. Note that for the antiproton-rich 2-gluon final state the improved and corrected 
formulas given in the second reference generally imply a lower branching ratio than those of the former reference which 
has been used in several previous analyses. 

The hadronization for all final states (including gluons) is simulated with the well-known particle physics Lund Monte 
Carlo program Pythia 6.115 (Sjostrand 1994), which is used extensively at accelerators in simulations of jet production 
at the full energy range which we need to consider here. A word of caution should be raised, however, that antiproton 
data is not very abundant, in particular not at the lowest antiproton lab energies which tend to dominate our signal. 
Therefore an uncertainty in normalization, probably of the order of a factor 2, cannot be excluded at least in the low 
energy region. 

3.2. Simulations 

To get the energy distribution of antiprotons for each of the final states listed in the previous subsection we generate 
the final states cc, bb, tt, W'^W~, Z^Z'^ and gg and let them decay/hadronize according to Pythia 6.115. We do not 
need to include lighter quarks since the branching ratios to these are negligible. The annihilation channels containing 
Higgs bosons need not be simulated separately since they decay to the other particles for which we do simulate. They 
are then let to decay in flight and the spectrum from the decay products are boosted and averaged over the decay angles. 
During the simulations, antineutrons are let to decay since we would otherwise underestimate the flux by a factor of two. 

We have performed simulations for the neutralino masses m^ = 10, 25, 50, 80.3, 91.2, 100, 150, 176, 200, 250, 350, 
500, 750, 1000, 1500, 2000, 3000 and 5000 GeV and for intermediate masses an interpolation is used. For each mass and 
annihilation channel, 2.5 x 10^ events have been simulated. 

For easier use, we have also parameterized the antiproton distributions for the 'basic' annihilation channels given above. 
It is for this purpose more convenient to deflne x = Tp/m^, with Tp being the kinetic energy of the antiproton, as the 
independent variable. A suitable parameterization is then given by 

^ = (pixP'^+P2iiogio^r)~' (4) 

where the parameters pi depend on both annihilation channel and the neutralino mass. The latter dependence is param- 
eterized as 

Pi{m^) = (a,iTO^'" + o,3m°*'')" (5) 

The values of the a^ for the different annihilation channels are given in Table 2. These parameterizations are valid for 
neutralino masses in the range 50-5000 GeV. The error in the most relevant regions, i.e. the low-energy tail to most of 
the high-energy slope is usually less than 20%. At worst (for tt), it can be up to 50% in isolated regions. This should be 



compared with the uncertainties in Pythia which probably caji be up to a factor of 2. In Fig. 2 we give as an example 
the distributions for a 100 GeV neutralino annnihilating into bb and a 1000 GeV neutralino annihilating into W'^W~ . 

Note that the above parameterizations are only given for the reader's convenience - in our calculations we use the results 
of the simulations directly. For convenience we also show in Table 3 the individual branching ratios of the nrain modes 
for a set of models with high antiproton yield but different mass, and gaugino content. 

3.3. The Antiproton Source Function 

The source function Q^ gives the number of antiprotons per unit time, energy and volume element produced in anni- 
hilation of neutralinos locally in space. It is given by 



)|(r,x) 



(fannl') 



Px(^) 



E 



dNf 
~dT' 



Bf 



(6) 



where T is the p kinetic energy. For a given annihilation channel /, B^ and dN^ / dT are, respectively, the branching ratio 
and the fragmentation function, and (cTannw) is the annihilation rate at f = (which is very good approximation since 
the velocity of the neutralinos in the halo is so low). As dark matter neutralinos annihilate in pairs, the source function 
is proportional to the square of the neutralino number density n^ = p-^/m^. Assuming that most of the dark matter 
in the Galaxy is made up of neutralinos and that these are smoothly distributed in the halo, one can directly relate the 
neutralino number density to the dark matter density profile in the galactic halo p. Given a generic parametrization of p, 
we fix: 



Px(^) = P(^) = Po 



ro 



7 r 



1 



(ro/a)" 



l + (|f|/a)" 



(/3-7)/q 



(7) 



where po is the value of the local halo density, ro is the galactocentric distance of the Sun and a is some length scale; we 
assume po = 0.3 GeVcm"^ and ro = 8.5 kpc. In the actual computation we will mainly restrict ourselves to the case in 
which the dark matter density profile is described by a modified isothermal distribution, (a,/3,7) = (2,2,0), mentioning 
what changes are expected in case more cuspy profiles, which are favoured by results in iV-body simulations of hierarchical 
clustering, are considered. We will in particular consider the example of the Navarro et al. profile (Navarro, Frenk & 
White 1996), (a,/?, 7) = (1,2, 1). Although we are here focusing on the case of a smooth distribution of dark matter 
particles in the halo, an extension to a clumpy distribution is potentially interesting as well (Bergstrom, Edsjo, Gondolo 
& UUio 1999; UUio 1999). 

Given the p distributions calculated in the previous section, we can now get the source function for any given annihilation 
channel. 
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Fig. 2. — Examples of the fits to the antiproton distributions from neutralino annihilation, dN/dT = {\ / m^)dN / dx with 
dN/dx given by Eqs. (4)-5) and Table 2. 



Example No. 


1 


2 


3 


4 


5 


6 


7 


M 


926.6 


-191.2 


708.7 


-127.6 


-543.1 


-891.2 


379.6 


Ma 


92.4 


-149.8 


-2156.1 


242.3 


413.3 


111.5 


-162.4 


tan/3 


19.1 


45.1 


9.67 


5.09 


4.17 


38.6 


28.8 


TUA 


79.8 


543.3 


1405.0 


1838.5 


409.2 


77.9 


89.5 


mo 


922.5 


306.6 


2817.0 


4243.2 


1054.6 


1110.6 


1109.1 


Ab/mo 


0.95 


0.17 


-2.05 


1.22 


-1.28 


0.17 


1.84 


At/mo 


-2.11 


0.96 


0.46 


-0.38 


-1.72 


2.06 


-2.20 


m^o 


46.2 


69.4 


704.9 


98.9 


208.2 


56.1 


80.9 


z. 


0.99761 


0.88843 


0.00642 


0.37943 


0.99293 


0.99756 


0.98440 


m^+ 


91.0 


117.9 


708.3 


123.8 


405.8 


110.9 


155.7 


rrijjo 


126.6 


543.3 


1405.2 


1838.9 


412.0 


126.8 


127.4 


rrijjo 


76.7 


107.5 


117.8 


115.5 


115.5 


76.2 


89.0 


m^± 


111.5 


549.7 


1407.1 


1840.1 


416.4 


115.5 


120.0 


(cTannW) 


0.127 


0.0874 


0.420 


0.0719 


0.116 


0.0824 


0.0785 


Br{cc) 


< 10--' 


0.00001 


<io-^ 


0.00050 


0.00007 


< 10-^ 


<io-^ 


Br{bb) 


0.96 


0.98 


0.81 


0.0059 


0.32 


0.96 


0.96 


Br{tt) 


- 


- 


0.11 


— 


0.65 


— 


- 


Br{W+W-) 


- 


- 


0.024 


0.87 


< 10"^ 


- 


< 10-^5 


BriZ^Z") 


— 


— 


0.020 


0.12 


< lO"'' 


— 


— 


Br{gg) 


0.00060 


0.0028 


0.00017 


0.0016 


0.0028 


0.00054 


0.00043 


Bri-yy) 


3.7 X 10"^ 


3.9 X 10"^ 


6.5 X 10"^ 


1.8 X 10"^ 


2.3 X 10"*^ 


9.6 X 10-^ 


1.2 X 10-" 


Br{Zy) 


5.8 X 10"^^ 


5.8 X 10"^ 


1.7 X lO""* 


6.5 X lO""* 


2.3 X 10"^ 


6.3 X 10-'° 


2.1 X 10-" 


n-^h'^ 


0.025 


0.028 


0.032 


0.025 


0.050 


0.034 


0.034 


$p 


4.6 X 10^2 


1.6 X 10"^ 


8.5 X lO-"* 


8.8 X 10-^ 


4.0 X 10"^ 


2.1 X 10-2 


1.1 X 10-2 


<&H, Earth 


3.8 X 10-' 


1.0 X 10° 


2.4 X lO"-"' 


1.1 X 10"^ 


2.5 X 10"^ 


7.5 X 10' 


2.3 X 10' 


<E>M, Sun 


7.0 X 10^2 


2.6 X 10° 


4.1 X 10"' 


2.7 X 10^ 


8.7 X 10-' 


1.3 X 10° 


1.8 X 10' 


O-SI 


2.9 X lO-'^ 


4.0 X lO"'^ 


9.6 X 10"'° 


1.2 X 10^* 


1.2 X 10-'° 


1.2 X 10-" 


2.3 X 10-" 


*e + 


5.0 X 10-* 


2.9 X 10"* 


3.2 X 10"* 


2.6 X 10"* 


2.0 X 10-* 


3.1 X 10-* 


2.6 X 10-* 


$cont. 7 


8.1 X 10"* 


3.5 X 10-» 


5.9 X 10-9 


1.4 X 10"* 


1.1 X 10-* 


4.2 X 10-* 


2.6 X 10-* 



Table 3 

Example of models giving high antiproton rates. All masses are given in GeV; the annihilation rate is given in IQ-^^ cm'^ 
s""'^; the solar modulated p flux at 0.35 GeV is in m"^ s~^ sr~' sr~^; the neutrino flux from the Earth and the Sun is 
with a threshold of 25 GeV and in units of km-^ yr-^; the spin-independent cross section is in pb; the solar modulated 



positron flux is in one of the HEAT bins (8.9-14.8 GeV) and in units of cm s 



GeV ; the continuum gamma flux 



is for high Galactic latitudes and is given in cm s sr , integrated above 1 GeV. The rates and fluxes are calculated 
as given in (Bergstrom & Gondolo 1996; Bergstrom, Edsjo, Gondolo & Uflio 1999; Baltz & Edsjo 1999; Bergstrom, Uflio 
& Buckley 1998; Bergstrom, Edsjo & Gondolo 1998). 



4. PROPAGATION MODEL 

In the absence of a well established theory to describe the interactions of charged particles with the magnetic field 
of the Galaxy and the interstellar medium, the propagation of cosmic rays has generally been treated by postulating a 
semiempirical model and fitting the necessary set of unknown parameters to available data. A common approach is to 
use a diffusion approximation defined by a transport equation and an appropriate choice of boundary conditions (see e.g. 
Berezinskii & al. 1990; Gaisser 1990 and references therein). 

In most diffusion models the form of the terms present in the diffusion equation is a compromise between physical 
insight and the possibility of an analytical solution. Only recently more realistic models have been studied by applying 
numerical solutions (Strong & Moskalenko 1998) or Monte Carlo simulations (Porter & Protheroe 1997). 

Our analysis is focused on comparing the characteristics of cosmic-ray antiproton signals of different origin: secondary 
antiprotons produced in cosmic-ray interactions and, eventually, a primary flux from neutralino annihilations. We want 
in particular to examine the dependence of the relative strength and spectral signatures on the diffusion model and the 
distribution of particle dark matter in the galactic halo. As both of them are not well constrained, we believe that 
an analytic solution of a reasonable physical model will be sufficient to provide most of the information needed on the 
behaviour of the two types of signals. 

We choose to describe the propagation of cosmic rays in the Galaxy by a transport equation of the diffusion type as 
written by Ginzburg and Syrovatskii (1964) (see also Berezinskii & al. 1990; Gaisser 1990). In the case of a stationary 
solution, the number density A'^ of a stable cosmic ray species whose distribution of sources is defined by the function of 
energy and space Q{E, x), is given by: 



dN{E,: 

dT 



= V • {D{R, x) \7N{E, f)) - V • {u{x) N{E, x)) - p{E, x) N{E, x) + Q{E, x) 



(8) 



Here and below we try to keep the notation as general as possible. Although our goal is to compute N for antiprotons, 
in order to determine the source function for the secondary flux we will have to obtain the spatial density distribution 
for protons as well. On the right hand side of Eq. (8) the first term implements the diffusion approximation for a given 
diffusion coefficient D, generally assumed to be a function of rigidity R, while the second term describes a large-scale 
convective motion of velocity u. The third term is added to take into account losses of cosmic rays due to to collisions 
with the interstellar matter. It is a very good approximation to include in this term only the interactions with interstellar 
hydrogen (on the other hand, we will point out below that heavier elements, in particular helium, cannot be neglected 
when computing the source function for secondary antiprotons) ; in this case p is given by: 

p{E,x)^n"{x)v{E)<jil^{E) (9) 

where n^ is the hydrogen number density in the Galaxy, v is the velocity of the cosmic ray particle considered ^cr\ while 
cr™ is the inelastic cross section for cr-proton collisions. In Eq. (8) we have neglected continuous energy losses; this will 
be included in an implicit form when considering secondary antiprotons. We will briefly mention in the conclusions the 
possibility that antiprotons have a finite lifetime r. To take this effect into account the term — (l/r) N{E,x) should be 
added on the right hand side of Eq. (8), and this corresponds to shifting p to p + l/r in all the equations below. 

We now have to specify the parameters introduced and the boundary conditions. We mainly follow the approach of 
Ginzburg, Khazan & Ptuskin (1980), given also in Berezinskii & al. 1990 and analogous to that of Webber, Lee & Gupta 
1992, Chardonnet & al. 1996, and Bottino & al. 1998. The main feature is that the propagation region is assumed to 
have a cylindrical symmetry: the Galaxy is split into two parts, a disk of radius Rh and height 2 ■ hg, where most of the 
interstellar gas is confined, and a halo of height 2 • hh and the same radius. We assume that the diffusion coefficient is 
isotropic with possibly two different values in the disk and in the halo, reflecting the fact that in the disk there may be a 
larger random component of the magnetic fields. We then have a spatial dependence: 

D{x) = D{z) = Dge{hg-\z\) + Dh9i\z\-hg) . (10) 

Regarding the rigidity dependence, fits to cosmic ray data in models which do not include reacceleration effects indicate 
that D scales as i?°'^ (Webber, Lee & Gupta 1992; Strong & Moskalenko 1998) with a cutoff below some rigidity Rq. We 
consider the same functional form as in Chardonnet & al. 1996 and Bottino & al. 1998: 



D,iR) = D\'^l + —j (11) 

where I ~ g, h. We will brieffy discuss below what changes are expected in case reacceleration is included, without making 
numerical predictions. The convective term has been introduced in Eq. (8) to describe the effect of particle motion against 
the wind of cosmic rays leaving the disk. We will therefore not consider any convection in the radial direction, assuming 
instead a galactic wind of velocity 

u(f) = (0,0,u(z)) (12) 

where 

u{z) ^s[gn{z)uhO{\z\-hg) . (13) 

An analytic solution is possible also in the case of a linearly increasing wind (UUio 1999). The last parameter we have to 
specify is the distribution of gas in the Galaxy: for convenience we assume that this has the very simple z dependence 

n^ix) = n"{z) - nf e{hg - \z\) + nf 6{\z\ - hg) (14) 

where Uh «C Ug (in practice we will take n^ = 0) and an average in the radial direction is performed. Finally, as we 
eventually want to treat the case of clumpy neutralino dark matter (UUio 1999) we will not assume any symmetry in the 
source function Q{x) (note that to apply the results of this Section to sources with a cylindrical symmetry, as for instance 
Eq. (6) with the assumption in Eq. (7), it is sufficient to set everywhere the index k equal to 0). 

As boundary condition, it is usually assumed that cosmic rays can escape freely at the border of the propagation region, 
i.e. 

N{RH,z)^N{r,hh)^N{r,-hH)^Q (15) 

as the density of cosmic rays is assumed to be negligibly small in the intergalactic space. To check whether this hypothesis 
holds even in the case of a source from dark matter annihilations, we will compare the flux of outgoing antiprotons with 
the one entering the diffusion region due to sources in 'free space'. For references about other possible choices of boundary 
conditions see Berezinskii & al. 1990. 

The cylindrical symmetry and the free escape at the boundaries makes it possible to solve the transport equation 
expanding the number density distribution iV in a Fourier-Bessel series: 

OO OO , \ 

n{t, ^, ^) = E E ^m ^']r ) ■ K'^^^ ^°'^^^^ + ^'^^^ ''""^H (^^^ 

k=0 s=l ^ ' 



which automatically satisfies the boundary condition at r = Rh, v^ being the s-th zero of Jk (the Bcssel function of the 
first kind and of order k). In the same way the source function can be expanded as: 

°° °° / r \ r n 

Q(r,z,0)-^^ J, /.,*—• Ql{z)cos{ke)+Q^,{z)sin{ke) (17) 

fe=0 s=l \ ^hj ^ 

where 

Rh -K 

Q'i^) = p 2 / 2, fc^ / d^' ^'Jk (ly^^) — I de' cos{k9') Q{r', z, 9') . (18) 

Rh Jfc+i {Vs) J \ RhJ akn J 

-7T 

In the equation above ag — 2, while ak — I for A: > 1; it is not necessary to specify the coefficients of the terms in shi(k9) 
as we fix the coordinate system such that 6* = at our location and we are only interested in computing fluxes for this 
value of 9. Inserting the two expansions in Eq. (8), we can derive the equation relevant for the propagation in the z 
direction: 

±Diz)-^^M^iz)-Diz)(^^^ M^iz)^^Ju{z)M^iz))-piz)M^{z) + Q';iz)^0 . (19) 

The solution of this equation is straightforward: it can be easily derived writing a solution separately for hg < z < hh, 
—hg < z < hg and —hh < z < —hg, and then imposing the boundary conditions at z = ±hh, i.e. Eq. (15), and the 
continuity of the number density and of the flux, that is of M^{z) and [—D(z)d/dz M^(z) + u(z) M^{z)], at z = ^hg. 
For —hg<z<hg the solution is given by: 

z 

M^{z) = Mi=(0)cosh(A,^^z) -^^J d^'sinh(Af (z - z')) Q^(z') (20) 

" '^ 



where 



^i^(O) - 7Zr7^^\-rTTlM-^r^ + ^^ b. + X'h'coth{X^,Hhh-hg))]+Iac 



COShiX^g^hg) \ Sinh (Af (/ift - hg)) Dg Af 



X [Dg \f tanh [Xfhg) +Dhlh + Dh A^^ coth {X^^hh - hg))] ' (21) 

and we have introduced the following set of definitions: 






and 



ks^, .,^^ „..„,., ,. .,^^ Q^(^') + 0^(-^') 



hh 

Ih = dz' sinh {X'^^'ihh - z')) exp {-fh{hg - z')) 

hg 

hg 

Igs = J dz' sinh {X'g%hg-z'))-9M±Mzll 





iGC 



dz' cosh (Af {hg -z'))- ^'^ ' '^'' '- . (23) 



For \z\ > hg the solution is analogous, but less compact and will not be reproduced here. It is more useful to give explicitly 
M^{0) in case of a source which is constant in the disk and negligible in the halo; in this case 



rkra, _ 1 It Dh [ih + ^t' COth {X'f,^ (hh - hg))] 



'^ ' Dg Xf^ \ Dg A^« sinh (Af /.g) + Dh cosh {Xfhg) [7. + Xf coth (A^^(/i, - hg))] J ^^^ ^ 



= M*(s,fc)Qj(0). (24) 

It is easy to check that in the limit Uh -^ 0, Eq. (24) correctly reduces to the result quoted in Berezinskii & al. 1990, 
Chapter 3, Section 3. Also Eq. (21) corresponds to the number density found in Chardonnet & al. 1996 and Bottino & 
al. 1998 in the limits hg -^ 0, Dh ~> Dg, Uh —^ and for a source function symmetric with respect to the z axis. 



Before applying these results to compute cosmic ray antiproton fluxes, let us pause for a moment. Most commonly, 
data on cosmic rays have been treated within the framework of the leaky box approximation. This is, to a certain extent, 
a simplified version of the diffusion model, where it is assumed that diffusion takes place rapidly. As was noticed for 
instance in Berezinskii & al. 1990, the path length distribution function of particles for a diffusion model with sources in 
the disk is very close to the exponential form characteristic of the leaky box treatment. For the purpose of computing 
secondary antiprotons, which are mainly generated in the gaseous disk, we expect it to be essentially equivalent to write 
a diffusion equation and fit its parameters to existing data on cosmic ray nuclei, or to derive from the same data in the 
simple leaky box scenario the grammage as a function of rigidity and use it to compute the secondary antiproton flux. 

In this sense we do expect to find an antiproton secondary spectrum which is analogous to the results of several papers 
in which this has been calculated in the leaky box approximation (except that in some of these papers not all the relevant 
effects have been included) . On the other hand, in the case of neutralino sources which are more homogeneously distributed 
extending through the full halo, it is unlikely that the effective average matter density particles have gone through can 
be of the same form as for sources located only in the disk. We will try to analyze in some detail this dependence on the 
geometry of the source, and we believe that this will give a real improvement with respect to the leaky box approximation. 

5. SOLAR MODULATION 

A further complication when comparing predictions of a theoretical model with data on cosmic rays taken at Earth is 
given by the solar modulation effect. During their propagation from the interstellar medium through the solar system, 
charged particles are affected by the solar wind and tend to lose energy. The net result of the modulation is a shift in 
energy between the interstellar spectrum and the spectrum at the Earth and a substantial depletion of particles with 
non-relativistic energies. 

The simplest way to describe the phenomenon is the analytical force-field approximation by Gleeson & Axford (1967; 
1968) for a spherically symmetric model. The prescription of this effective treatment is that, given an interstellar flux at 
the heliospheric boundary, d^b/dTb, the flux at the Earth is related to this by 



pI dn 



^(^-)-3^(^^) (^^) 



where the energy at the heliospheric boundary is given by 

Ei, = Es> + \Ze\<j)F (26) 

and Pig, and pb are the momenta at the Earth and the heliospheric boundary respectively. Here e is the absolute value of 
the electron charge and Z the particle charge in units of e (e.g. Z = —1 for antiprotons). 

An alternative approach is to solve numerically the propagation equation of the spherically symmetric model (Fisk 
1971): the solar modulation parameter one has to introduce with this method roughly corresponds to (pp as given above. 
When computing solar modulated antiproton fluxes, the two treatments seem not to be completely equivalent in the low 
energy regime (the reader may check for instance Fig. 4 in (Labrador & Mewaldt 1997) against Fig. 8 in (Bottino & al. 
1998)); keeping this in mind, we will anyway implement the force field approximation, avoiding the problem of having to 
solve a partial differential equation for each of our supersymmetric models. 

We just mention here that the picture can be much more complicated: non-spherical propagation models which take 
into account the polarity of the solar magnetic field have been studied as well (e.g. Webber & Potgieter 1989). In this 
case the solar modulation effects on particles and antiparticles can be quite different. If one would translate this into the 
force field effective treatment, one should use different values for the modulation parameter for protons and antiprotons 
(in contrast with the standard procedure of assigning to antiprotons the value found from proton flux measurements). 
The relation between these two would be very model dependent. 

To compare with the two sets of Bess measurements, which are both near solar minimum, we choose (pp = 500 MV, in 
reasonable agreement with what the Bess collaboration uses in their analysis. Otherwise, we will focus on the interstellar 
fluxes which are not affected by these uncertainties. 

6. BACKGROUND ESTIMATES 

6.1. General considerations 

Secondary antiprotons are produced in cosmic ray collisions with the interstellar gas. Looking at the composition of 
incident and target particles it is easy to guess that the main contribution to the p flux is given by cosmic ray protons 
colliding with interstellar hydrogen atoms. Because of baryon number conservation, the minimal p + p ^f p -\- X reaction 
has three protons in the final state; in the rest frame of the target hydrogen atom, the energy threshold for an incident 
proton to produce an antiproton is therefore Ep =^ Imp. Due to this feature, the energy distribution of the produced p 
shows a sharp peak at a few GeV and a steep fall-off at lower energies. It is reasonable to expect that in this low energy 
region reactions involving heavier nuclei, both as targets and projectiles, may play some role: they imply in fact different 
kinematics and the spectrum of the produced p need not fall as fast as for low energy pp collisions. We have verified 
that the interaction of primary protons with interstellar helium is indeed a relevant process, while all others can be safely 
neglected as their contributions add up to below a few per cent of the total at any energy (this is essentially the same 
conclusion as that of Simon, Molnar & Roesler 1998, although our approach is slightly different, as we will point out 
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below). We assume therefore that the source function for secondary antiprotons has the foUowing form (the factor of 2 
accounts for antiprotons produced by antineutron decays): 



CO 

]p{x, S) = 2 • 47r f dE' 

-^thresh 



^^I^iE,E')n"ix) + ^^^^{E,E')n"%x) 



%ix,E'). (27) 



In this formula ^p{x,E') is the primary proton cosmic ray flux at the position x in the Galaxy and for the energy E', 
n^'^ is the helium number density which we assume to be 7% of n^ (Garcia-Munoz & al. 1987) and have the same spatial 
dependence, while da/dE {E, E') stands for the differential cross section for the production of an antiproton with energy 
E for an incident proton of energy E' , in the two processes considered. For pH collisions we implement the standard 
parametrization for the differential cross section introduced in Tan & Ng 1983; as already mentioned the energy threshold 
-E'throsh for this process is E' = Imp. 

6.2. Antiproton production in collisions with nuclei 

A much smaller set of data is available in the case of antiproton production in proton collisions with heavier elements. In 
particular, it is diflicult to estimate the effects of production below the nominal energy threshold for p + p ^ p + X, which 
is known to occur in hadron-nucleus collisions. Recently, several experiments have shown substantial sub-threshold p 
production for deuterium, helium, carbon and copper targets (Chiba & al. 1993; Schroter & al. 1993). Possible collective 
eflfects allow the abundant low-energy cosmic rays to produce antiprotons below threshold. On their way out of the 
nucleus, the produced antiprotons may also suffer inelastic losses which slow them down, creating a potentially important 
component in the low-energy cosmic ray-induced p spectrum. (Note that helium and heavier nuclei in the cosmic rays 
also give different kinematics for produced antiprotons. However, this gives an extra contribution at higher p energies, 
and is therefore not important for our study.) These sub-threshold effects have been modeled by Sibirtsev & al. 1997, 
where the limited data set available can be described by a transport equation solved by a Monte Carlo technique. 

Here we follow a much simplified approach, which describes the C and Cu data displayed in Sibirtsev & al. 1997 
reasonably well and which we apply to p + He ^ p + X. We find that the collective effects can be mimicked by a shift 
in the incident proton energy: 

Ein -^ Ecff — -Ein + 0.6(ii'throsh — -Ein)^(£^throsh — -^in) + 1-1 GcV, (28) 

where i^throsh is the nominal threshold {Imp) for p production in pp collisions. The energy loss due to inelastic rescattering 
can be approximated by decreasing the energy of the outgoing antiproton (using pp kinematics) by 1 GeV. The yield of 
antiprotons per collision is taken to scale with the total pA cross section, parametrised according to Letaw, Silberberg & 
Tsao 1983. In this way, we have a parametrisation which is asymptotically correct at high energies and which also fits the 
subthreshold data. However, we are unable to assess the accuracy of this treatment for the problem at hand, believing it 
to be at the 50 % level, but acknowledging the need for improved data and theoretical modeling of this seldom discussed 
problem. 

6.3. Primary proton flux 

The last step to make before implementing Eq. (27) is to determine the primary proton flux '^p{x, E'). The formalism 
introduced in Section 4 is suitable for this purpose once we specify the source function for primary protons. It is generally 
believed that supernova remnants are the main sources of cosmic rays. Nevertheless the gradient of $p as a function of 
the distance from the galactic centre which is obtained from the observed distribution of supernovae or the related pulsar 
distribution is not consistent with models for gamma-ray emission (see Strong & Moskalenko 1998, and references therein). 
We take advantage of the phenomenological approach of Strong & Moskalenko 1998, where a generic form for the radial 
distribution of cosmic-ray sources was considered and its parameters fitted to EGRET gamma-ray data (Eq. (6) and 
Fig. 12 in Strong & Moskalenko 1998). We therefore assume that the primary proton source, in cylindrical coordinates, 
is of the form: 

Qp{E,x)^q{E)q{x)^q{E)(^\ C^^ (-'—I^\ e{h, - \z\) (29) 

where rg = 8.5 kpc is our galactocentric distance, and we have assumed that the energy spectrum of emitted protons 
is the same everywhere in the Galaxy. The function q{E), which we may interpret as a normalization factor, can be 
rewritten, after propagation, in terms of the local proton flux $p(ro, E) which has been measured in several experiments. 
It was argued in the past that the spread among different experimental determinations of <I>p (rg , E) introduces one of the 
main factors of uncertainty in the prediction for the secondary antiproton flux (see for instance Gaisser & Schaefer 1992). 
The recent measurements by the Imax (Menu & al. 1997) and Caprice (Boezio & al. 1999) collaborations are in better 
agreement with each other. In Bottino & al. 1998, a fit of the data of these two experiments with a single power law in 
energy or rigidity was made, using the force-field method to take solar modulation into account. The fits in rigidity (see 
also Boezio & al. 1999) show a steeper fall-off than those in energy; however, this may not be the case if a break in the 
spectrum at low rigidities is assumed instead (Ormes & Protheroe 1983). As in the next Section the background at high 
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energies will be important, we prefer to conservatively consider fits with a power law in energy. From Eq. (1) and Table 1 
in Bottino & al. 1998, 

M^o,i^)^AvL^^^__j (30) 

with A = 12300 ±3000 and a = 2.67 ±0.06 for the Imax data, and A ^ 19600 ±3000 and a = 2.85 ±0.04 for the Caprice 
data. 

6.4. Interstellar secondary antiproton flux 

We are now ready to give the formula for the interstellar secondary antiproton flux at our galactocentric distance. We 
find: 

<^p{ro,E) = ^vp{E)Np{ro,E) = 



9 fi7^^ T f o^o\ ,,. p. , . . f ,^, ^{E,E')n^%{ro,E')Mj,{s,E') 
= ^'"p{E)}^Ju[v,—-]Mp{s,E)lR{s) / dE -^ ^ 

where the repeated index i stands for the sum over hydrogen and helium, and we have introduced the notation: 

Rh 



(31) 



Rh'Ji'{^^)J \'Rh 



Ir(^) = TT^irn^ I dr' r'Jo ( ^."^ ) q{r') (32) 







and Mp{s,E) — M*(s,k — 0) setting the total inelastic cross section appropriate for pp collisions ct^" = cr'" , while 
Mp{s, E) — M*{s, k — 0) with (t™ = cr!f . The parametrizations for both of these cross sections were given in Tan & Ng 



1982: 



<p(r)=<( 32.2- 
32.2- 



32.2.[l+0.02731n(g/200)] , "^ < T ^ -^ PpV 

ri+o.oo262-T-("8+i3**i''T+4-4ii„aT)| mt) , u.d S J < ^ t^ev 



1 + 0.0273 ln(^/200)] mb , 3 < T < 200 GeV (33) 

1 + 0.0273 ln(£;/200) + 0.01 (ln(.B/200))^ 



mb , T > 200 GeV 



and 



(7fp{T) = 24.7 • (1 + 0.584 T-o"'"^ + 0.856 T^O'^^*^) mb , T > 0.05 GeV (34) 

where T and E are in units of GeV. To derive Eq. (31) we have assumed that for z < hg the approximation ^p{z) c± 
$p(z = 0) is valid. This is generally a very good approximation, as for most choices of the parameters in the propagation 
model (f>p is nearly constant in the disk and rapidly decreasing in the halo (see for instance Fig. 3.10 in Berezinskii & al. 
1990). Only in extreme cases can ^p{z — hg) be 10% lower than $p(z = 0) and the correction to the result in Eq. (31), 
always below few per cent, can be obtained keeping track of the full z dependence in Jp (use Eqs. (20) and (21), all 
numerical integrals in Eq. (23) can still be performed analytically; the result follows easily). 

In Bottino & al. 1998, it was suggested that it is a good approximation to assume that the energy spectrum for the 
protons after propagation is roughly independent of location in the Galaxy. This hypothesis simplifies the computation 
(for us, it is especially needed to compute numerically the tertiary contribution described below); Eq. (31) reduces to: 

OO 

2vp{E) J dE'^{E,E')n^%{r^,E') ^ 
'^pinuE) = ^^^ ^ ^ Y. Jo ( ^°^ ) Mp{s, E) Mp{s, E) Ir{s) (35) 



E7=i J,[^^Jt)Mp{s',E)In{s') 

where E is an arbitrary normalization energy. As there are some indications that the energy spectrum may indeed be 
steeper far away from the sources, because of the energy dependence in the propagation coefficient (Mori 1997), we 
compare in one case Eq. (35) against Eq. (31) to check if the simplification in any way changes the result. For the set of 
parameters as in example in Fig. 1, which we will soon discuss, we find that Eq. (35) gives a very slight overestimate of 
Eq. (31), below 3% for interstellar antiproton kinetic energies up to 1 GeV, a maximal 5.5% overestimate at 3 GeV, while 
for higher energies the difference decreases again and is below 4% at 50 GeV (we remark however that we have not tuned 
our propagation model to reproduce the effect in Mori 1997 so we cannot claim that this effect is not relevant). 

6.5. Tertiary antiprotons 

In Eq. (8) we have not introduced any energy-changing term. We now include energy losses for secondary antiprotons 
due to scattering processes during their propagation in the Galaxy. The main effect is due to non-annihilation inelastic 
interactions of antiprotons with interstellar protons, giving lower energy antiprotons in the final state. Actually, the 
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Fig. 3. — The interstellar antipron flux and the contribution from secondary and tertiary antiprotons. The uncertainty 
due to the parametrization of the primary proton spectrum is also given as the shaded band. The solid line corresponds 
to the same set of parameters as in Fig. 1. 

energy distribution after the non-annihilation interaction is not well known as there are no direct measurements; the usual 
assumption (Tan & Ng 1982) is that the distribution is similar to the final state proton in pp inelastic (non-diffractive) 
interactions, i.e., a rather flat distribution in kinetic energy between zero and the kinetic energy of the incident antiproton. 
One may think that elastic scattering processes are relevant as well but available data show that the cross section is 
dominated by the forward peak with very small energy transfer (Eisenhandler & al. 1976; Briickner & al. 1986), and hence 
with a marginal net effect from our point of view. We include therefore only non-annihilation processes considering a 
'tertiary' source function generated by inelastically scattered secondary antiprotons in the form: 



Qf\x,E) = ATrn''{x) 



noil— ann 

Upp 



{E') 



T' 



Ip{x,E')dE' 



non — anil 
(Tpp 



{E)Ip{x,E) 



(36) 



where a^™ '*"" is obtained as the difference between the total inelastic cross section Eq. (34) and the inelastic annihilation 
cross section: 



pp 



'(r) = 



661 • (1 + 0.0115 T-o-^'^^ _ 0.948 T00151) 
36 r-"-5 mb 



nib 



T < 15.5 GeV 
T > 15.5 GeV 



(37) 



where for lower energies we have used the parametrization in Tan & Ng 1982, while in the high energy range we apply 
the approximation given in Protheroe 1981. Both this parametrization and those needed above have been checked against 
a compilation of more recent data (Caso & al. 1998). The second term in Eq. (36) takes into account antiprotons which 
are depleted from the energy E and which we propagate as a negative flux; it actually gives an effect that is less than 
few per cent at any energy and is not needed in our formalism. As was done for Eq. (27), it is straightforward to write a 
Fourier-Bessel expansion for Q^"''* and then compute ^p""^* which has to be summed to $p to get the final expression for 
the background interstellar antiproton flux. 

6.6. Numerical results 

Coming to the actual numerical predictions for the background flux of antiprotons, we base our choice of parameters 
in the propagation model on previous work in which diffusion models analogous to the one described in Section 4 were 
used to fit data on cosmic-ray nuclei, such as ratios of secondaries to primaries and of radioactive nuclei to their stable 
counterparts. Actually slightly discrepant results are present in the literature, partially refiecting the fact that it is not 
easy to find a propagation model which is consistent with the whole set of existing data. We consider here and in the 
following Section, when describing the signal from neutralino annihilations, three different scenarios. We only keep the 
following parameters fixed: n^ = 1 cm~'^, nj^ — 0, hg ^ 0.1 kpc and Rh = 20kpc. The first three are the standard values 
which are inferred from direct observation. The last one, which in the literature is taken sometimes as small as 15 kpc. 
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Fig. 4. — The effects of changing the diflFusion parameters are shown. In (a) the thin halo scenario, in (b) the thick halo 
scenario, while in (c) a convective halo is considered. See the text for further details. 



and which in Strong & Moskalenko (1998) is set equal to 30 kpc, does not play a major role and different choices lead to 
nearly equivalent results. 
The three scenarios are: 



a) 



In the case D" 



Dl 



D^ our propagation model is fairly close to the one considered by Webber, Lee and Gupta 
(1992). Their conclusion is that a thin halo is preferred, with height hh G (1-1, 3.8) kpc and D ~ (6±4)-10^^ crn^s"^ at 
the rigidity R= 1 GV. The flux shown in Fig. 1 is obtained in this scenario, setting hh = 3 kpc, Z?° = 6 • 10^^ cm^s~^ 
and i?o = 3 GV, choosing the proton flux at the Earth as the medium value in the fit of IMAX data, i.e. with 
A = 12300 and a = 2.67, and taking into account solar modulation with the force fleld method with (f>p = 500 MV 
as suggested by the analysis of the Bess collaboration. There is no consensus in the literature on the value of the 
solar modulation parameter at solar minimum, nevertheless the spectrum does not change dramatically if a slightly 
different value for (j)p is assumed. For instance, ipp — 400 MV gives about a 7% increase at the kinetic energy 
T ~ 0.2 GeV and about an 8% increase at the maximum. For (pp = 600 MV the effect is reversed and wc find that 
the flux is lower by roughly the same percentages as in the previous case. 

As will become clear in the following, the background antiproton flux shown in Fig. 1 is only an example of the 
possibility of a good flt to the data; we keep it as reference case to compare with. In Fig. 3, we show for the 
same parameters the interstellar antiproton flux versus kinetic energy T, plotting also its three main components: 
the secondary antiproton flux due to pp collisions, the contribution from pHe scattering processes and the tertiary 
component due to energy loss. As can be seen the first contribution is dominant at the maximum and at high 
energies, while the other two are important in the low energy region. 

We take advantage of Fig. 3 to show another feature that is common for all choices of the propagation parameters, 
the uncertainty due to the interstellar proton flux. The band around our reference antiproton flux is the envelope of 
the predictions obtained by using the uncertainty in the proton flux (Bottino & al. 1998): the upper bound is given 
choosing the fit of IMAX data with A — 15300 and a = 2.61 (average values +1ct and —la respectively), while the 
lower bound below T — 2.5 GeV is obtained from the IMAX data fit with A = 9300 and a = 2.73 and above 2.5 GeV 
from the CAPRICE data fit with A = 16300 and a = 2.89 (average values —la and +la respectively, actually these 
two spectra are nearly overlapping at all energies). 

Coming back to the uncertainty in the choice of the propagation parameters in the Webber-Lee-Gupta scenario, if 
we now pick the average value for the halo height h^ = 2 kpc and vary the diffusion coefficient in the suggested 
interval, I?° ~ (3 — 7) • 10^^ cm^s~^ for Rq = 3 GV, we find that the the flux at intermediate energies increases by up 
to about 30% for the smallest value of the diffusion coefficient, while it decreases with a slightly higher percentage 
for the highest value of D°. This is represented by the band in Fig. 4, region (a) (both in this Figure and in Fig. 6 
below we defined fractional differences as ($ — $/{)/$/{, with $/j being the reference value, i.e. in this case the flux 
shown as a solid line in Fig. 3). If we assume on the other hand that D" and hh are linearly related, a degeneracy 
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which may indeed not be resolved by available data, and fix D'^ = 2.5(/ift/kpc) ■ lO^^cm^s"^ varying hh between 
1.1 and 3.8 kpc we get a band of very small width, below a few percent (solid and dashed lines at about —20% in 
Fig. 4, region (a)). 

b) In Strong & Moskalenko 1998, a scenario is favoured with a thicker halo, with hi^ G (4, 12) kpc, in case no convection 
is assumed. Their treatment of propagation is less close to our model than the previous case, so the way we translate 
their typical choice of parameters into our picture is less safe but should give at least the right qualitative behaviour. 
We sketch the thick halo scenario taking hh G (4, 12) kpc, Df^ — 2.5(/i/i/kpc) • 10^"^ cm^s^^, D^ ~ 6- 10^^ cm^s^^ and 
i?o = 1 GV, a choice consistent with the results in Ginzburg, Khazan & Ptuskin (1980), where the propagation model 
we have chosen was first considered. Comparing to our reference case we find a band of 20% width around an average 
suppression of the flux of about 50%, where the least severe suppression is given by the smallest halo considered and 
the maximal suppression corresponds to the large halo h^ — 12 kpc (solid and dashed lines respectively in Fig. 4 
region (b)). 

c) As a third scenario we allow for the presence of a galactic wind driving cosmic rays out of the galactic disk. Self 
consistent models for the propagation of cosmic rays in magnetohydrodynamic flows have been studied recently 
in detail (Zirakashvili & al. 1996; Ptuskin & al. 1997). The much simpler approach we take here is intended to 
compare qualitatively the effects of the wind on the background antiproton flux and on the signal from neutralino 
annihilations. Considering again the model in the previous scenario with /i^ = 4 kpc (solid curve in Fig. 4 region 
(c)), we take as an example the case of Vh = 10 km s^^ (lower dashed line) and Vh — 20 km s~^ (lower dash-dotted 
line). As can be seen, in perfect analogy with the case of the solar wind, the galactic wind alters the spectrum at low 
energies (up to 30% in the example we are considering) while the effect gets smaller and smaller for more energetic 
particles. If in analogy with the parameter choice of Strong and Moskalenko we suppose that there is some scaling 
between w/i and D^, for instance a simple linear scaling D° = 2.5(/ift,/kpc)[(40km s~^— Vh)/40km s^^] -lO^^cm^ s~^, 
we find that at intermediate energies the flux is nearly unchanged (upper dashed and dash-dotted lines). 

We will not combine the uncertainty bands we have just derived and make a deflnite statement about the uncertainty 
on the prediction of the cosmic antiproton background. To be able to do that on a firmer basis we should compare the 
predictions of our propagation model directly against the whole set of data on cosmic ray nuclei and this is beyond the 
aim of the paper. We stress again that this Section was mainly intended to show that the most recent data on cosmic ray 
antiprotons can be fitted by the background flux for some natural choice of the diffusion parameters. On the other hand 
we find that the prediction for the background could be lower as well, leaving room for an antiproton flux generated by 
an exotic source, possibly dark matter ncutralinos. 

6.7. Reacceleration 

It seems very plausible that cosmic rays are reaccelerated by a Fermi type of acceleration by stochastic magnetic 
fields during propagation. This has been treated, e.g., in Seo & Ptuskin (1994); Heinbach & Simon (1995) and Simon 
& Heinbach (1996). There is also a possibility that cosmic rays get reaccelerated by weak shock waves from supernova 
remnants (Letaw, Silberberg & Tsao 1993). We will here focus on the former process, usually called diffusive reacceleration 
since it can be treated as a diffusion in momentum space. 

In Heinbach & Simon (1995) it was shown that data on low energy cosmic rays are compatible with the predictions of 
models without diffusive reacceleration only if the mean path length variation with energy shows a sharp break around 
1-2 GeV. In models which include reaccelaretion effects on the other hand, depending on reacceleration strength, a path 
length distribution that is a simple power law for all energies may be considered. This is theoretically appealing since this 
from observations derived form agrees well with that expected from Kolmogorov turbulence. The result in Heinbach & 
Simon (1995) is confirmed in the analysis by Strong & Moskalenko (1998), who conclude as well that the reacceleration 
scheme allows a more natural choice for the parameters in the propagation model. 

Even though reacceleration implies that the average energy increases as the cosmic rays propagate through the galaxy, 
there is a general smearing of the injected spectrum meaning that there is also a 'leakage' of cosmic rays from higher 
energies to lower. 

This might be important for antiprotons where the injected spectrum drops below the maximum at few GeV. In Simon 
& Heinbach (1996) it was found that this leakage could substantially increase the secondary p spectrum below 1 GeV, 
with the flux at a few hundred MeV that hardly can be lower than 1/3 of the value at around 1 GeV. We cannot compare 
directly with our analysis as we are not applying the same primary proton spectrum and propagation model, but still 
we can conclude that including reacceleration in the propagation model might add up to the effects of the proton-nuclei 
interactions and of the tertiary component in flattening the antiproton spectrum at low energies, and making it even more 
problematic to separate an exotic signal from the background in this region of the antiproton spectrum. 
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Fig. 5. — The value of Cprop for the same choice of diffusion parameters as in Fig. 1 and an isothermal sphere distribution 
of dark matter with a — 3.5 kpc (solid line), compared to the value for the same diffusion parameters and a Navarro et al. 
profile with a ~ 9 kpc (dotted line) . The band gives the range of values of Cprop in the thin halo scenario (with h^ = 2 
kpc) as described in the text. 



7. SIGNAL FROM NEUTRALINO ANNIHILATION 

7.1. General discussion 

With the source function introduced in Eq. (6), the antiproton flux from neutralino annihilations in the galactic halo is 
readily obtained from the formulas derived in Section 4. It is given by: 



47r 47r ^ — ' 

s— 



>o(.."^)mO(0) 



(38) 



where M°(0) is obtained from Eq. (21) with crj.^p = a^^^ in Eq. (22) and Q°(z) in Eq. (23) given by 



Rh 



(39) 



It is possible to separate in the expression for the signal the part which depends on the MSSM parameter space from 
the terms which are related only to the distribution of sources in the propagation region and to the propagation model 
itself. We introduce the definition: 



dNf „f / po ^ ^ 



<^piro,T) EE (a,n„z;) ^ -—Bf ^ Cp,op(r) 



(40) 



The quantity Cpiop, which can be obtained explicitly by comparing Eq. (38) with Eq. (40), has the dimension of length 
divided by solid angle and is analogous to the coefficient defined in Eq. (46) of Bottino & al. 1998; note however that in 
Eq. (40) we have factorized the value of the local halo density po rather than some reference density. 

7.2. Uncertainties related to propagation 

Having factorized out in Eq. (40) the dependence of the signal on the choice of the dark matter candidate, we analyse 
first how sensitive the result is to the set of parameters which define both the location of the sources and the propagation of 
the produced antiprotons. We fix a reference configuration selecting for the propagation model, in analogy to the analysis 
of the background fiux, the same parameters as in the example in Fig. 1 and Fig. 3 {h/^ = 3 kpc, Z?" = 6 ■ lO^^cm^s""'^ 
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Fig. 6. — The changes of Cprop when the diffusion zone height is changed within (a) the thin halo scenario derived from 
Webber, Lee and Gupta model and (b) the thick halo scenario inspired by the Moskalenko and Strong analysis. In (c) 
the effects of a convective wind are sketched: in (cl) the diffusion coefficient is kept fixed while in (c2) it is linearly scaled 
with the galactic wind. 



and i?o = 3GV), while as a reference dark matter density profile we choose a modified isothermal distribution, Eq. (7) 
with (a,/3, 7) = (2,2,0), and with an intermediate value for the length scale, a = 3.5 kpc. In Fig. 5 we plot the value of 
C'prop for this reference case (solid curve) versus the antiproton kinetic energy. Cprop is increasing in the low energy range 
as it contains the kinematic factor vp{T) and at the same time we have assumed that the diffusion coefficient is roughly 
constant at low rigidities (see Eq. (11)). The coefficient Cprop then reaches a maximum at about T ~ 2 GeV, while at 
higher energies it decreases as a consequence of the 0.6 power law increase in the diffusion coefficient. 

First we analyse how the result changes if as the dark matter density distribution we consider the Navarro et al. profile 
(Navarro, Frenk & White 1996) which is singular towards the galactic centre, Eq. (7) with (a,/3,7) = (1,2,1). Choosing 
a = 9 kpc, we obtain (see the dotted line in Fig. 5) roughly a 33% increase in Cprop (and therefore in the signal from 
neutralino annihilations) at any T, while the more cuspy profile with a = 3.5 kpc gives a result which is more than 
twice the reference value. This is rather surprising because even though the singularity in the profile induces a sharp 
enhancement in the neutralino number density and therefore in the strength of the source, this cusp is at the galactic 
centre rather far away from the solar system. It is commonly believed that in the diffusion regime the local sources are 
the most relevant, but at least in the propagation model we are considering, this is not true: for the Navarro et al. profile 
with a = 9 kpc, 23% of the signal is given by sources contained in a spherical region of 1 kpc around the galactic centre; 
the percentage increases to 42% for a = 3.5 kpc while it is as low as 1% for the isothermal sphere profile. We conclude 
that non-local sources may give a significant contribution provided their strength is much enhanced with respect to the 
local ones. This effect should be considered in more detail when considering a clumpy scenario (UUio 1999) for which it 
may be even more relevant. 

Fig. 5 contains another piece of information. Going back to the case where the dark matter density profile is described 
by an isothermal sphere, we have varied the parameters which define the propagation model, going through the same three 
scenarios described when discussing the background. The band in the figure is given by fixing h^ — 2 kpc and varying the 
diffusion coefficient in the interval D° ~ (3 — 7) • lO^^cm^s"^, and corresponds to the band shown in part (a) of Fig. 4 
(again the highest value of D'^ gives the lowest value for the fiux). Unlike the latter, the band in Fig. 5 does not overlap 
the reference value (solid line): for the signal from neutralino annihilations the decrease in the height of the propagation 
zone from 3 kpc to 2 kpc is not compensated by the decrease in the central value for the diffusion coefficient. This gives 
a first hint on how sensitive the dark matter signal is to the choice of the value of the height of the diffusion zone. 

The same effect is studied in the upper part of Fig. 6, fixing the kinetic energy to T = 1 GeV, varying hh and linearly 
relating D^ to hh, as introduced in the previous Section. According to the Webber-Lee-Gupta scenario, hh is constrained 
to be between 1.1 and 3.8 kpc: the degeneracy we found for the background flux (the nearly overlapping solid and dashed 
line in part (a) of Fig. 4) is completely removed for the signal from neutralino annihilations, going from an 80% suppression 
to a 20% increase compared to the reference value, for minimal and maximal h^ respectively. Changing the diffusion zone 
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Fig. 7. — In (a) the interstellar antiproton flux at 1 GeV is shown and in (b) the ratio of the flux at different kinetic 
energies to that at 1 GeV is shown. To make the flgure clearer (and avoid showing artifacts of sampling frequency) the 
figure is binned. We also show with different symbols in which bins there are models which are gaugino-like, mixed and 
Higgsino-like. 



height modifies the number of sources which contribute to the flux, as sources which arc outside the diffusion box are not 
included in the model. Therefore it is not surprising that a very thin diffusion zone gives a suppressed signal while larger 
values for hh enhance it. 

In the same way, in the thick halo scenario without convection, inferred from the analysis of Strong and Moskalenko, 
the suppression band found for the background flux (part (b) of Fig. 4) becomes a much wider band for which the signal 
flux is increased instead (solid line in part (b) of Fig. 6). The enhancement with the diffusion zone height is flattened 
out at high values of h^ as the new sources we include are further and further away from the observer and moreover the 
density proflle falls at large galactocentric distances. 

While the effects we have considered so far give roughly the same result for any value of the antiproton kinetic energy, 
the effect of convection in the z direction is clearly energy dependent. In part (cl) of Fig. 6 we plot Cprop as a function 
of the galactic wind speed, for hh = 4kpc, D^ — lO^^cm^s"^, Z?° — 6 • lO^^cm^s^^ and Rq — 1 GV (we consider here 
as reference value the one obtained with this set of parameters and Uh ~ 0). In part (c2), as we did for the background, 
D^ is taken to be related linearly to the value of the wind speed. In both parts flve different kinetic energies (given in 
the Figure in GeV) have been considered. As can be seen going back to part (c) of Fig. 4, the effect of convection is 
greater on the signal from neutralino annihilation than on the background flux. Especially at low energies, it is not well 
compensated by the linear scaling of the diffusion coefficient. 

For the background flux the effect for introducing a galactic wind is to drive antiprotons more quickly from the disk, 
where they are generated, to the border of the diffusion zone where they are lost. This effect can be balanced by lowering 
the diffusion coefficient, that is, assuming that diffusion takes place less efficiently. The sources of the signal are on the 
other hand distributed over the whole diffusion box; setting with the galactic wind a preferred direction of propagation 
lowers the probability for an antiproton generated relatively far away from the disk to reach our location. The effect is 
not compensated by the rescaling of DjJ, at least not for the rescaling needed for cosmic ray species generated in the disk. 

The last check we perform regards the role played by the antiprotons which are produced outside the propagation 
region. Our solution to the diffusion equation has been derived under the hypothesis that the number density of the 
considered cosmic ray species is zero at the boundary of the diffusion zone. This is not strictly true for the signal from 
neutralino annihilations. One possibility of verifying what kind of corrections might be needed in this case is to compare 
the antiproton flux leaving the diffusion zone with the flux injected by external sources. We restrict the analysis to 
exchanges at the boundary z = i/i/i, as the effect is much suppressed in the radial direction, being generally Rh » hh- 

The outgoing flux can be computed by keeping track of the full dependence on z of the number density, as this flux 
is related to the gradient of the number density at the boundaries; one derives a rather lengthy expression which we do 
not reproduce here but which follows in a straightforward way. For the injected flux we use the very simple picture of 
propagation in free space, summing contributions over the line of sight. For the modified isothermal sphere profile and a 
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Fig. 8. — The interstellar antiproton flux at 1 GeV versus the relic density. This is the only figure where models having 



fi^/i^ < 0.025 are shown. 



diffusion zone height of 3 kpc we find that the ingoing fiux is about one third of the outgoing fiux for very small radial 
coordinates, while they become roughly equal at r = 8 kpc and at larger radii the injected fiux becomes prevailing; the 
total number of antiprotons per second that penetrate the diffusion box is about 70% of those that leave it. 

As can be understood, this fraction is smaller if we consider instead the Navarro et al. profile or if we pick a higher 
value for the diffusion zone height, and might turn into a rather large number for very thin haloes. To give a precise 
numerical estimate for the effect one should add in the propagation model a third zone above \z\ — hh- However, taking 
into account all the other uncertainties that enter in the prediction for the signal, we do not consider this worthwhile at 
present. We believe that a safe assumption is that the signal from neutralino annihilations has not been underestimated 
due to this effect by more that a factor of 2 in the most extreme cases. 



7.3. Antiprotons from specific MSSM models 

The antiproton spectrum from neutralino annihilation have been calculated for all the different MSSM models given in 
Table 1. In Figs. 7-11 we show our main results. We use our canonical parameters for propagation and the isothermal 
sphere model for the halo profile. 

7.3.1. Interstellar fluxes 

In Fig. 7 (a) we show the predicted interstellar antiproton flux at a kinetic energy of 1 GeV (i.e., without corrections 
for solar modulation) versus the neutralino mass. We clearly see the trend that the fiux goes down with the mass of the 
neutralino. The reason for this is that the number density of neutralinos goes down as the mass increases, for a given 
interval of the dark matter mass density. Since n^ = p^/m^, and the annihilation rate scales as n^ the suppression 
increases rapidly with mass. At the lower mass end, the present accelerator limits preclude a neutralino in the MSSM 
below a few tens of GeV. The low-flux models at low masses will be discussed in connection to Fig. 8 below. 

The points in the figure are coded with different symbols for different composition of the neutralino. We define models 
with Q < Zg < 0.01 as Higgsino-like, 0.01 < Zg < 0.99 as mixed and 0.99 < .Z^ < 1 as gaugino-hke. As can be seen, most 
of the models with high rates are either gaugino-like or mixed, except at masses greater than several hundred GeV, where 
also Higgsino-like models can be important. (In fact, there is also a small mass window around 80 GeV where Higgsinos 
may be relevant.) 

In Fig. 7 (b), we show the ratio of the interstellar flux at 0.5, 3, 5, and 10 GeV to the fiux at 1 GeV displayed in (a) for 
the same set of models (but without coding the composition). It is seen that for the higher mass models, it can be more 
advantageous to study the fiux at higher kinetic energies. 

In Fig. 8 we show the same fiuxes as in Fig. 7 (a), but versus the relic density, Vt^h^ . There is a very clear trend that 
the highest flux is obtained when Q.^h? is close to the lowest acceptable relic density. The reason for this is that if the 
annihilation cross section is increased, the fiux of antiprotons increases, but since the relic density Q,^h^ is approximately 
inversely proportional to the annihilation cross section, Q,^h? decreases and hence the strong correlation. The correlation 
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Fig. 9. — The flux of antiprotons from neutralino annihilation at the optimal kinetic energy, Topt, versus Topt- Topt is 
defined as the energy at which $signai/'&background is highest and if the spectrum has more than one optimum, the highest 
two have been included in the plot. The models have been coded according to the neutralino mass in GeV. 



is not perfect however, since it is the thermally averaged cross section at a temperature of about TOp^/20 that determines 
rt-^h^, whereas the annihilation in the halo to a very good approximation occurs at rest (the speeds are typically ~ 0.001c). 

This also exlplains some features in Fig. 7 (a). In the mass range between 40 and 60 GeV there exist models which 
give exceedingly small rates, but also some which give the highest of all rates. This large spread reflects peculiarities 
near the Z° (and neutral Higgs) resonances and the W^W~ threshold. For the low-flux models around 40-60 GeV, the 
annihilation cross section at rest is very small, but either a resonance or threshold can be reached through thermal motion 
in the early Universe and the relic density is reduced to our selected range 0.025 < Vt^h"^ < 1. As shown by Chen & 
Kamionkowski (1998), three-body final states can be important (not too far) below the W^W~ and the ii thresholds and 
this could enhance the signal for these low-rate models. 

We also have some models around 130 GeV that give high fluxes. In this case it is the other way around, the masses 
are just so that we are on the iJg-resonance for the non-relativistic speeds in the halo, but the thermal average in 
the early Universe gives a lower annihilation cross section and hence the relic density is increased to our desired range 
0.025 < n^h^ < 1. The behaviour at 130 GeV is just accidental, it could happen at any important resonance. In fact, we 
found only one high-flux model around 130 GeV in our 'normal' scans, and performed a small scan varying the parameters 
slightly around this model. The relic density was essentially unchanged, but the antiproton flux showed large variations 
depending on if we were below, on or above the resonance. 

In Fig. 8 we also show models with a value of flh^ lower than our required limit of 0.025. In principle, one could accept 
these models at the expense of introducing other components of dark matter. To be consistent, one should then rescale 
the local dark matter density in the form of neutralinos by some unknown factor. In the lack of better procedures, one 
usually employs a linear rescaling p^ = (11^/1^/0.025) -pDM- Since the annihilation rate is quadratic in the number density, 
this rescaling factor enters squared in the predicted p rate, something which is clearly visible in Fig. 8. 

We are now interested in finding out if there are any special features of the antiproton spectra from neutralino annihi- 
lation which distinguish these spectra from the background. We have already mentioned that the window at low energies 
may not be as good as previously thought. We will here investigate other features and energy regions of the spectrum to 
see if there are good signatures of a neutralino contribution to the flux. 

One thing that might differ is the slope at different energies. The background is expected to have a rising trend at low 
energies, reaching a maximum between 1 and 2 GeV (see Fig. 3) and a slope of around —3 at high energies. On the other 
hand, the high-rate models tend to be decreasing at 1 GeV (see Fig. 11 (a)). This may cause a shift of the maximum of 
the summed spectrum (signal plus background) to a lower energy, which is a possible signature. 

We next investigate if there is an optimal energy at which $signai/$background has a maximum (for this purpose we will 
use the reference background given in Fig. 3). In Fig. 9 we show the flux at these optimal energies, Topt, versus Topt. We 
now see that we have two classes of models. One class which have highest signal to noise below 0.5 GeV (i.e. inaccessible 
in the solar system due to the solar modulation) and one which have highest signal to noise at 10-30 GeV. For this first 
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Fig. 10. — The solar modulated antiproton fluxes at (a) 0.35 GeV and (b) 1.76 GeV compared with Bess 97. The models 
have been coded according to their relic density, Q^h"^. 



class of models, we note that there exist a proposal of an extra-solar space probe (Wells, Moiseev & Ormes 1998) which 
would avoid the solar modulation problem and is thus an attractive possibility for this field. However, these models have 
high rates in the range 0.5-1 GeV as well, even though it would be even more advantagous to go to lower energies. The 
second class of models are much less affected by solar modulation and also give reasonably high fluxes. In Fig. 11 (a) we 
show some examples of spectra. All of these have optimal energies in the low-energy region, but e.g. spectrum 3 has an 
optimum at high energy as well. 

7.3.2. Solar modulated fluxes 

We now turn to the solar modulated fluxes and will compare with the Bess 97 measurements, which we recall are in 
very good agreement with our estimate for the background flux. We will compare the fluxes in two of the Bess energy 
bins, the one at 0.35 GeV and the one where the measured flux is the highest. In Fig. 10 we show the solar modulated 
fluxes versus the neutralino mass. We see the same general trend as for the interstellar fluxes. Fig. 7, but we also see that 
there are many models with fluxes above the Bess measurements. However, this conclusion depends strongly on which 
range one allows for the neutralino relic density. In Fig. 10 we have coded the symbols according to the relic density 
interval. As can be seen, essentially all models which are in the Bess measurement band have a relic density fl^h^ < 0.1. 
If we instead require 0.1 < Q^h^ < 0.2 the rates are never higher than the measured flux. 

This points to a weakness of this indirect method of detecting supersymmetric dark matter: once the predicted rate is 
lower than the presently measured flux, the sensitivity to an exotic component is lost. This is because of the lack of a 
distinct signature which could differentiate between the signal and the background. Alternative indirect search methods, 
like gamma rays from the halo (see e.g. Bergstrom, Ullio & Buckley 1998), or neutrinos from the Sun or the Earth (see 
e.g. Bergstrom, Edsjo & Gondolo 1998) have the added virtue of giving both a directional and a spectral signature which 
can be used to improve the signal to background ratio well beyond the limits of present-day measurements. 

The highest values for the fluxes in Fig. 10 (a) are 4 times higher than the Bess measurement. However, the uncertainty 
coming from the local halo density alone is larger than this. Given the total mass of the galaxy, and restricting to our 
choice of halo profile (isothermal sphere with a — 3.5 kpc), we find a minimal local halo density of 0.14 GeV/cm^ which 
would correspond to a flux reduction of a factor of 4.6. To that one should add the uncertainties of the Monte Carlo 
simulations (up to a factor of 2) and the halo profile, the propagation model and solar modulation. For this reason it is 
presently not possible to exclude any supersymmetric model on the basis of antiproton measurements alone. 

7.3.3. Example of models 

In Table 3 we show 7 MSSM models that all give high p fluxes. These models have acceptable relic densities, cover a 
large mass range and have varying composition (and obey present accelerator bounds). 

In Fig. 11 (a) the predicted differential p flux is shown for the 7 models. They show maxima occurring at lower energies 
than for our canonical background. At higher energies, the trend is that the slope of the flux decreases as the neutralino 
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Fig. 11. — (a) Antiproton spectra for all 7 models appearing in Table 3. (b) Example of a composite spectrum consisting 
of our reference background p flux (Fig. 1) reduced by 24 % with the addition of the predicted flux from annihilating dark 
matter neutralinos of MSSM model number 5 in Table 3. 



mass increases. Model number 3 corresponds to a heavy neutralino and its spectrum is significantly less steep than the 
background. If such a spectrum is enhanced, for instance by changing the dark matter density distribution, we would get 
a bump in the spectrum above 10 GeV. 

In Table 3 we also show the annihilation rate and the most important branching ratios. With the help of the results 
in the earlier Sections, the antiproton flux from neutralino annihilation can be derived. The difference between the 
parameterizations given in Section 3.2 and the full simulation results (that we have used) is typically less than 20%. 

Also note that the branching ratio to gg is never important for our high-rate models (not only the ones in the table, 
but all high-rate models). This is not in agreement with the results found by Jungman & Kamionkowski (1994). The 
reason for this difference is that we use the improved gg annihilation cross section of Bergstrom and Ullio (1997). 

The table also contains an indication of the rates for other detection methods. The neutrino-induced muon flux in 
neutrino telescopes does not show a strong correlation with the p flux, and it is possible to find models that give either low 
or high rates. Current limits are about lO'^-lO^ muons km~^ yr~^. We also give the spin-indepedent neutralino-nucleon 
cross section (Bergstrom & Gondolo 1996), which should be compared with the current limits that are of the order of 
10~^ pb. These show a better, but not perfect, correlation with the p fluxes. The correlation is even stronger between the 
p flux and both the e+ fiux and the 7 flux with continuum energy spectrum. Both of these do not decrease as much with 
neutralino mass as the antiproton flux does, however. For more details, see Baltz & Edsjo (1999) and Bergstrom, Edsjo, 
Gondolo & Ullio 1999. The cross section for annihilation into monochromatic 7s (through 77 and Z7) are uncorrelated 
with the p flux. 

Finally, in Fig. 11 (b) we show an example of a hypothetical composite spectrum which consists of our canonical 
background flux decreased by 24 % (obtained e.g. by decreasing the primary proton flux by la), and the signal for model 
5 in Table 3. We can obtain a nice fit to the Bess data, but as noted before, there are no special features in the spectrum 
that allow us to distinguish between this case and the case of no signal. 

8. DISCUSSION AND CONCLUSIONS 

We have seen that there is room, but no need, for a signal in the measured antiproton fluxes. We have also seen that 
the optimal energy to look for when searching for antiprotons is either below the solar modulation cut-off or at higher 
energies than currently measured. However, there are no special spectral features in the signal spectra compared to the 
background, unless the signal is enhanced and one looks at higher energies (above 10 GeV). 

We have stressed the somewhat disappointing fact that since the present measurements by the Bess collaboration 
already exclude a much higher p flux at low energies than what is predicted through standard cosmic-ray production 
processes, an exotic signal could be drowned in this background. Even if it is not, the similar shape of signal and 
background spectra will make it extremely hard to claim an exotic detection even with a precision measurement, given 
the large uncertainties in the predicted background flux (at least a factor of a few, up to ten in a conservative approach) . 
We note that some of the uncertainties may be reduced if an extrasolar probe aimed at low-energy detection would be 
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launched, a possibility that has been recently proposed (Wells, Moiseev & Ormes 1998). 

Although it is tempting to conclude that what has been measured by the Bess experiment is the standard cosmic-ray 
induced background flux of antiprotons, one should keep in mind that it could, on the contrary, be almost entirely due 
to an exotic source like neutralino annihilation. Since this possibility cannot be excluded (at least until the problem of 
the dark matter in the Galactic halo has been solved), one has to be cautious about using the measured antiproton flux 
to deduce properties of antiproton propagation and, as has recently been done (Geer & Kennedy 1998), the antiproton 
lifetime. We have checked that, using one of our high-mass neutralino models and a clumpy distribution of dark matter 
in the halo, we can get an excellent fit to the Bess data for antiproton lifetimes as low as 10^ years, clearly violating the 
claimed lower bound of Geer & Kennedy (1998). (For details, see UUio (1999).) 
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